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Zeros of the Riemann zeta function and the resulting distribution of the primes 
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Abstract: The Riemann zeta function is one of the most Euler’s important and fascinating functions in mathematics. 

By analyzing the material of Riemann’s conjecture, we divide our analysis in the C(z) function and in the proof of the 
conjecture, which has very important consequences on the distribution of prime numbers. The proof of the Hypothesis 
of Riemann result from the simple logic, that when two properties are associated, (the resulting equations that based in 
two Functional equations of Riemann ), if zero these equations, ie C(z) = £ (l-z)= 0 and simultaneously they to have the 
proved property 1-1 of the Riemann function C(z).Thus, there is not margin for to non exist the Re (z) = 1/2 {because 
C (z) = C (l-z)=0 and also tj(z) as and 1/(1-z) are 1-1}.This, as it stands, will gives the direction of all the non-trivial 
roots to be all in on the critical line, with a value in the real axis equal 1/2. 


Part I, Proof Riemann Hypothesis 


y l — y 2 = 2kzr /ln(ft), is an integer multiple of In. But 


#7. Theorem 1. 

The R-Hypothesis focuses on the point where we 
must prove that if s = Re(s) + Im(s)*I... 

I) The functions ^(s) and £(l-s) are 1-1 on the critical 
strip. 

II) The common roots of the equations £(s)-£(l-s) = 0 
they have Re(s)=Re(l-s)=l/2 within the interval (0,1) 
and determine unique position, which is called 
critical line. 

Proof: 

The functions £(s) and £(l-s) are 1-1 on the critical strip. 
For this we need to analyze when and where the 
exponential function n z is 1-1 when neZ + . 

The exponential function n z ... 

I) is 1-1 in each of these strips defined by the intervals 
27ik/ln(n) and 27i(k+l)/ln(n) “ where k e Z . For this we 
must prove two cases.. .If AcCa/:AaC 

then if /(z 1 ) = /(z 2 ) then = Z 2 that is.. 
x 1 = x 2 a y 1 = y 2 • Indeed, if z x = x x + i • y x , 

£ 2 = x 2 + i • y 2 are two points within into a such strip 

such that e Zl =e Zl then,, 

=> | n ‘i| = |„*a| => ^ n y^ = y 2n y2^ => 

X 1 = X 2 and since X 1 = X 2 , the relation 


Zi and Z 2 belong to the strip, so \y t - y 2 \ < 2zr/]n(n) 

That is, the difference y - y 2 is at the same time a 
multiple of 2nAn(n) and at an absolute value of less than 
27i/ln(n).The only case that this is true is when 

y i = y 2 . We finally conclude that z 1 = z 2 , son z is 

1-1 in the strip 

{z : 2Jbr/ln(n) < Im(z) < 2(k + l)^/ln(n)} .’’lines 

down closed - open up”. We also notice that n z is on 

C-{0}.Because if w ^ 0 and we put 
z = ln I w\ /]n(w) + /argw/ln(w) then 

^z w \ /ln(n)+iargw/ln(n) ^In|w|^/arg w/\n(n) II ^iargw/ln(n) 


2^/Utfi 


2 (k + l)*ri/Ln(n) 


X\ + iyi 


■ Lvi - y*\ < /Ln(n) 


x 2 + % 
2km /Ln(n) 


Fig.l The Z 1 and Z 2 belong to the lane width 2 n / In (n) 


Formation of the strips 1-1 for ^(s) and ^(1-s). 
al.) If we accept the non-trivial zeroes on critical strip 
of the Riemann Zeta Function £(s) as s 1 =cr l -\- it l and 


Zi-ln(n) _ z 2 -ln(n) 


i-y i ln(^) i-y ? -ln(«) 

gives e 1 = e y2 , so 


s 2 = cr 0 + it 2 with \s 2 1 > \s { |, and if we suppose that the 
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real coordinates cr 1 , cr 2 of each non-trivial zero of the 
Riemann Zeta, [1,2,3,4] function £(s) correspond to two 
imaginary coordinates t\ and t 2 , then, we have the 
following equations group: 


C(cr1+1-^) = 


°i+^i 






- + ....+ - 


o x +vt x 


- + .... = 0 


C, (c2 + l' t 2 ) —- — -^--^- Z~t -- ^ --l~-“ 

~ v ^cr 2 +l-t2 2 <T 2 +i ' ? 2 ^ cr 2 +i ' ? 2 ^cr 2 +t-t 2 

Taking the first equation and deducting the second, we 
obtain: 


00 1 1 

£(°i+*-O- ~ ^ a 2 + 1 ■ h) = £(-^- -^^7) = 

n =1 n n 


00 M °2 + *‘ r 2 M Cr l + '' r l 00 

Z Vl —VI ^ 


— yf 2 e itrln(jl) — n a] e itr]n(n) 


= 0 


n=1 n 


From the previous relation, we conclude that if 

n a2+l ' t2 . n a i +l '' f i ^ 0, > 1 then 


2kn 

<J 2 -G i At 2 -t x ±-,(/; = 1,2,...) .That is to say ti 

ln(ft) 

and t 2 can take any value, but according to the previous 
relation. So from (Fig.l), that means £ (s) is 1-1 on the 

lane of critical Strip in Intervals such as defined, 

therefore and on the critical Line. So rf z is 1-1 in the 

strip and{z : 2kn!\ciri) < Im(z) < 2(k+ V)7r/\n(n)} 

.’’lines down closed - open up” 

a2.) If we do the same work with |s 2 1 > \ s \ | for the case 


£(l-s) will we have: 


00 1 1 


® n '-°2-‘h -rf-vi-'h 
L l-o- 2 -^2 

n =1 n n 


00 n~ a2 e ~ u ^ n) — e ~ u ^ n(n) 


n 


1 -<72-1*2 . 


= 0 


From the previous relation, we see that if 

n i-<r 2 -it 2 , n ^-‘-h ^ 0, n > 1 

&({l-cr 2 = 1 - cr,} v cr 2 = cr,) a 

_ 2kn /r . _ . 

At i =t x +-—-Xk = lL...) 

ln(» 

So, from (Fig.l), it follows that £ (1-s) is 1-1 on the 
critical Strip, So n” (1_z) is 1-1 in the strip and 


{z : 2k7r/]n(n) < Im(z) < 2{k + V)7r/]n(n)} .’’lines up 
closed - down open” 

As we have seen before, the complex exponential form 

n z but and n z is also periodic with a period of 

27i/ln(n). For example n z+2 ' K ' 1 = n z • n 27r ' l/Log{n) - n z the 
n z is repeated in all the ’’lines down closed - open up” 
horizontal strips with 27i/Ln(n) on (Fig.2) below. 



+6iTi/Log(n) 


+4TTi/Log(n) 


+2TTi/Log(n) 


+2TTi/Log(n) 

° - y 


-2ni/Log(n) 


-4Tri/Log(n) 


-6TTi/Log(n) 




Fig 2. All the "bottom - closed upper “horizontal strips 
with 27i/ln(n) for rf z , n > 1, n e Z . 


If we are asked to find the strips that are 1-1 of the 

function ^(^) = V— it will be the union of the strips 
7ln z 

formed with period 2 • n • k / In n where e Z . 

This results from the analysis we made for the cases al, 
a2. Therefore, if we analytically assume that n = 2, 3, ... 
then the exponential complex function will be 1-1 in 
each of the strips, i.e. in each one set of the form... 

S n = {z : 2k7r/\n(n) < Im(z) < 2(k + l)7r/\n.(ri)},k e Z 

and therefore in each one subset of the unified 

expanded set.. . S = Uj^ 

The largest strip obviously contains all the rest it is wide, 
and has width 2 • n / In 2 i.e. with n = 2. 

Overlapping contiguous Strips. 

If we start from the upper bound of the strip with n = 2 
and k = 1, Fig.3 the upper ones will descend to zero by 
increasing n with a width difference the interval 5 = 
27i/Ln(n-l) - 27i/Ln(n). In addition, we will prove that the 
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strip i=n>=2 primary (k=1) that generated in ascending 
order, the upper limit of that created of its double width 
enters the zone of the previous of (I = n-1, k = 1) also 
primary.. 











o 


Fig3. The upper ones will go down to zero 
by increasing n. 

Apply 2-7T<2-2-7T/]nn=> n<(n-l) 2 => 


ft 2 -3-ft + l>0=>?z> 


3 +V 5 


and because we talked 


about integers then n > 3 .The same phenomenon is also 
created at higher level strips. 


The 




i'=i » 


it is directly 1-1 from the 


comparison of 2 points 


Zl,Z 2 


, when 


C(z 1 ) = C(z 2 ) 


Proof: 

From the equality 

^) = I 4 = I-W(Z 2 ): 


_ _|_ 2-(U+yH) +yU) _|_ n ~(x\+y\-i) _ ^ _ 


_ Y _|_ 2~ ( - x2+y2i ) _|_ y(x2+y2i) j^-(x2+y2-i) _ s 

with comparison of similar terms the result arises. 


*1 = * 2 A y x = y 2 ■ 


2 '7i 'k 
In (n) 


,k e Z; 


If we need z 1 = Z 2 => = x 2 a y 1 = y 2 therefore it will 

be 1-1 as shown in sections al, a2 pages 1-2... Also 
according to Lagrange's generalized theory, each of the 
roots according to the obvious relation ^ (z) — C,{\ — z) 
will result from the generalized theorem of Lagrange. 
From the equation that results from 3 Functional 
Equations, with condition of common roots. We take 
4'(z) = £(l — z) and for first approach we then follow 


the 3 first terms on each ^-equation we will have in the 
Regular form £(z),£( 1 — z) with the analysis below. 

r + 2~ z + 3" z = r°- z) + 2^ + 3 H1 ~ Z) => 2" z + 3’ z = - 2 Z +- 3 Z => 

2 3 

2" z + 3 _z -2" <1_z) -3“ <1 - z) =0=>6 Z + —4 Z -3-3(-) z = 0 

2 3 

With the replacement 

t t , \ i • Log(y) + 2-k-7T-i 

6 =y=>z-Log(6) = Log(y) + 2-k-7r-i^>z = --——- 

Log( 6) 

_ Log(y)-h2-k ■ n -i 

Z Z ~ Log 6) 

we consider it as initial value for the solution 
of Transcendental equation with the method of 
Lagrange inversion, then: 


_ Log(y) + 2-k-7r-i 
Log (6) 


Hr 


jr(w + l) dy' 


T« 


--■)■(—■ 4 

y-Log{ 6) 2 


Log (6) 


-3-A 


Log(y)+2-k-n-\ 
Log{ 6) 


n, 


Some solutions for k £ Z (Table. 1) 

1st 0.4985184869351 +/—3.6072014300530 I,k = 1 
2st 0.4997077392995 +/-10.329023684711 I, k = 3 
3st 0.5003426560276 +/-13.990255078243 I,k = 4 
4st 0.4993087825728 +/-17.585326315022 I,k = 5 
5st 0.5026747758005 +/- 21.304192633721 1,k = 6 
Table.1-6 Roots of equation £(z) = £(l — Z ) 

And so on for the sequence of infinity. We now 
see that we approach the real roots and its complexities. 
With the additional method Newton or the Bisection 
method we approach more rapidly to them. And these 
roots will be contained in corresponding strips that we 
have previously defined and in closer relation with their 
position on the critical line. 

By this logic for the previous approximation equation in 
relation to its solutions of, ^(z)=0 the strips which 
determine the imaginary roots of the equation, will be 
between the intervals 2-;r-/f/ln6,/feZ as defined by 
the analysis of the generalized Theorem of Lagrange. 

The strips of the imaginary part of the roots, as they 
appear for 6 consecutive intervals of the approximate 
equation with method of Lagrange that we have solved 
above and we will have the results on the (Table.2), 
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Strips - Imaginary Roots: 


1st 
2st 
3 st 
4st 
5st 
6st 


(7.013424977055179 

(10.52013746558276 

(14.02684995411035 

(17.53356244263794 

(21.04027493116553 

(24.54698741969312 


3.506712488527589) 

7.013424977055179 

10.52013746558276) 

14.02684995411035) 

17.53356244263794) 

21.04027493116553) 


Table.2- The strips of the imaginary part of the roots 


It is clear the first root 3,607..i is located in 1 st of strip 
(7.01-3.506), the second 10.3290..i in the 2 st strip 
(7.013-10.52) etc. Now the root 13.99..i which 
approximates the first root of £(z) = ^(1-z) = 0 i.e. the 
14.1347..i, is in the 4 st strip (14.02-17.533) with a lower 
limit the 14.026..i a value that is very close to the one 
required for the approximate and exclusive (1 st root) when 
£(z) = 0 of the Riemann Hypothesis. 

Corollary 1. 

The only roots of the Zeta function not included in the 
set {z EC: 0 < Re (z) < 1} are the points -2, -4, -6.. 

[H.p 47] 

Proof.. 

On the functional equation... 

4X1 - z) = 2 • £(z) • r(z) • cos(/r/ 2 • z) • (2 n)~ z we know that 

for Re (z)> 1 the functions £(z) and T(z) do not equal 
zero(Proof-Th.3.p6).We also find that for Re (z)> 1 the 
1 -Re (z) <0 and putting u = 1 -z we will find all the roots 
of £ (u) for Re (u) <0.Therefore we will have that £ (u) 
will be zero , where it is zeroed the cos(;r/ 2 • z) i.e. 
z = 3,5,7..Then the roots of £ (u) for Re (u) <0 will be the 
points -2, -4, -6... and will be all the roots of the function 
4 out of the strip {z EC: 0 < Re (z) < 1}. 


Theorem Helpful. 

The Riemann’s Z-function has no roots on the lines 
Re(z) = 1 and Re(z) = 0. 

The proof is detailed in the book [11. p50-51]. 


II.)The common roots of the equations 4(s) - 4(1-s) = 0 
they have Re(s)= Re(l-s)=l/2 within the interval (0,1) 
if moreover apply 4(s)=0.(Refer to Theorem 3, p.5). 
Proof: Let us assume z to be such that for complex 

Z = x 0 + (y 0 ,O<Re(z) <lAlm(z) ^0 and 


4Xz) = (TO-“ Z) = 0. According to the two equations, 
they must apply to both, because they are equal to zero 

that... £( Z ) = 4Xz) = 0 a £(1 - z) = 4X1-7) = 0. But 

from [Theorem 1, II, al, a2] the 4(z) and 4 (1-z) are 1-lon 
the critical Strip.. If suppose generally that 

£(x 0 ± y 0 -i) = £(l-x' 0 ± y' 0 -i) = 0 then: 


To = To A *^o — 1 — x o ’ but because are the 1-1 then we 
will apply two cases for complex z i.e. 
x o =1- x o>x 0 ^Xo,y 0 =y' 0 v 


vXg =1 -x 0 ,x 0 =x' 0 =3,y 0 =y' 0 

[Corollary.l and Theorem{Helpful},because we 
suppose £(z) = (T(l — z) = 0. ] We conclude that 

0< X 0 <1. Therefore we have three cases: 


a) If 



we apply the obvious i.e. 


ai/2±y 0 -i) = ai-i/2±yo-i) = ai/2±yo-i) = OA 


a y 0 = Yq • Which fully meets the requirements of 


hypothesis!! 
b) 


o<x 0 *x;*-<i 


In this case 


if x 0 < Xq a x 0 + Xq = 1 or x 0 > Xq a x 0 + x' 0 =1 
for the Functional equation ^(s) - ^(1-s) = 0 
simultaneously apply: 


then 


^(X 0 +y 0 i ) = ?(X 0 -y 0 ' ) = 0 Aftx'o +y' 0 -i ) = 


= £(Xq — y c i) = 0 • So let's assume that... 

1 st .X 0 <Xg aX 0 +Xq =1=>x 0 =l-aAXo =^ + b , a * b . 

But apply x 0 +x[,=1=>1 + b- a = 1=>b = a. That is, 

a,b is symmetrical about of 1. But the Qz) and C, (1-z) 

2 

are 1-lon the critical Strip then apply 







5 


^-«±yo- i ) = 4'(^+«±yo-i)Ay 0 = yo}=> must 

— — a = — + a => 2a = 0=> a = 0 

2 2 


therefore 



2 nd . x 0 > Xq a x 0 + Xq = 1 => x 0 = ^ + a a Xq = ^ - a. 

as before. Then because the £(z) and £ (1-z) are 
1-1 on the critical Strip will apply... 

C(^+a±y 0 -i) = ^-a±y' 0 -\)Ay 0 =y / 0 } that 

—+ a = — —a => 2a = 0 => a = 0 
2 2 


therefore 



.We see that in all acceptable cases 


is true that “The non-trivial zeros of £ (z)- £ (l-z)=0 
have real part which is equal to Vi within the interval 
( 0 , 1 ) • In the end, we see still one partial case that 
appears in the roots of the equation 

az)-ca-z)=o. 


1 < X 0 A Xq < 0 


Furthermore, in this case withthat 


is, it is symmetrical about 


X Q + Xq = 1 then for the functional equations we apply 


£( x o _ yo •') = £(i- x o - y 0 ••) a £(x 0 +y 0 - i) = 


with regard to the three cases (a, b,c } the cases (b, c} 
cannot happen, because they will have to be within (0,1) 
and will therefore be rejected. 

Therefore if C(z) = C(l~z) = 0 then because u these 
two equations of zeta function are 1-1 on the lane of 
critical strip” as shown in theorem 1, is “sufficient 
condition ” that all non - trivial zeros are on the critical 
line Re(z)=l/2. That means that the real part ofz of 
C(z) = 0, equals to 1/2 So Theorem 1& Corollary 1. 
they have been proved and also helps Theorem 3 below. 


2#. Theorem 2. 

For the non- trivial zeroes of the Riemann Zeta 
Function £(s) apply 

i)There exists an upper-lower bound of Re(s) of the 
Riemann Zeta Function £(s) and more specifically in the 
In 2 In n 


closed interval 


In 2 n In 2n 


. The non-trivial zeroes 


of the Riemann Zeta function £(s) of the upper-lower 
bound are distributed symmetrically on the straight 
line Re(s) = 1/2. 


ii).The average value of the upper lower bound of 
Re(s) = 1/2. 


Proof: 

i). Here, we formulate two of the functional equations 
from E - q. Set 


=C(i- x o+yoi) A y 0 = Jq and for the three cases it is 
valid y 0 =Ovy 0 ^ 0) . This case applies only when 


Zfz. ) = £(1 - Z ) ^ 0 a X 0 = 1 + a, Xq = -a , and not for 

the equality of £ function with zero. But in order to 
have common roots of the two functions £ (z) and £ (1-z) 
«because they are equal to zero and equal to each 

other», then first of all, it must be i.e x 0 = Xq = 1 - X 0 

and also because the £ (z) is 1-1 and for each case 
applies (y Q = 0 and y' Q = 0) or (y 0 * 0 andy^-*= 0), 


1 

therefore we will apply Xo=Xo=-Ay 0 =yo 


and 


£ (1- s) /£ (s) = 2(2n )~ s Cos (ns / 2)r (s). Res >0 
C (s)/t; (l-s) = 2(2tt )" _1 5m(TO/2)r(l- sfRes <1 

We look at each one equation individually in order to 
identify the set of values that we want each time, 
a). For the first equation and for real values with 
Re(s) > 0 and by taking the logarithm of two sides of the 
equation [5], we have... 

C(1 -s)/£(s) = 2(2xy s Cos(7isl2)T(s) => 

Log[C(1 ~ s) / = Log[ 2 ] - sLog[2n\ + 

+ Log[Cos(nsl2)T(s)\ + 2 km but solving for s and if 
f(s) = Cos(7Ts/2)F(s) and from Lemma 2[6,p6] 

if Vm£(l-s)/£(s) = l or 

s^>s 0 

Log[ lim C( 1 - s)/£(s)\ = 0 we get: 

s^>s 0 
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Log[2] | Loglf(s)} + 2km with 
Log [27r] Log\l7i\ 

„ / L0g[/(5')] + 2^\^ „ r: i i i 

Re(-) > 0, jinaly because we need 

Log[ In] 

real s we will have Re(.y) > = o 3771 . This is 

Log[ 2n] 

the lower bound, which gives us the first £(s) of 
Riemann’s Zeta Function. 

b). For the second equation, for real values if Re(s) < 1 
and by taking the logarithm of the two parts of the 
equation, we will have: 

C(s) / C(1 -s) = 2(2 Try- 1 Sin(ns / 2)r(l -s)=> 

Logics) / ca - S)] = Log[ 2] + (s- l)Log[27T] + 

+ Log[Sin(7T-s/2)r(l-s)] + 2kmbut solving for s 

and if / (v) = Sin(7is I 2)T(1 — s ) and from 

Lemma 2[6,p.6] if lim £(s)/£(l-s) = 1 or 


This suggests for our absolute symmetry of 



\ o- 0 + it 7 

\ 

ai + iY 0 * 1 

1 •<x 2 +if 0 

N 

0 

Res = 

<7 0 + ttl 

=y 2 a 

» 


Fig.4-Arrangement of low and upper of real part of 
zeros from 2 functional equations of Riemann. 

ii) The average value of the upper lower bound is 
Re(s) = 1/2 because from (fig. 4): 

ReW=1/2 .if«ra±i£sw = i/2 

_ Log[27T\ _ 

#3. Theorem 3 

The Riemann Hypothesis states that all the nontrivial 
zeros ofC ( z ) have real part equal to 1/2. 


Log [lim C(s) / £(1 -s)] = 0 we get... 

Login] Log[f(s)] + 2km with 
Log[ In] LogYln] 

R U, g [f(s)] + Wm 
Log[ In] 

In the following, because we need real s we will take.. 

ReO) < = 0.622854 -This is the upper 

Logl 2n] _ 

bound, which gives us the second of Riemann’s Zeta 
Function of £(s). So, we see that the lower and the upper 
bound exist for Re s and they are well defined. 

c). Assuming that s k = cr low + it k and s ’ k = a upper + it k 


with (7 

lo\ 


L °g [2] and 
Log[2n] 


G 


upper 


LogW 

Log[2n] 


we 


apply cr, = <J low and cr 2 = a upper ,(FigA). If we 


evaluate the difference and 

ReAs; =Re.;-l/2 = a wer -l/2 = 0.1228 
ReAy K =l/2-Re5 K =l/2-a tow =0.1228 


Proof: In any case, we Assume that: The Constant 
Hypothesis £ (z) = 0, z ^ C. In this case, we use the 
two equations of the Riemann zeta function, so if they 
apply what they represent the £ (z) and £ (1-z) to 
equality. Before developing the method, we make the 
three following assumptions: 

I). Analysis of specific parts of transcendental 
equations, which are detailed... 

a) . For z e C ,{a z = 0^><2 = 0a Rez > 0 } which 
refers to the inherent function similar to two Riemann 

zeta functions as (2 7l)~ z — 0 or (2 7 z) z ~ l = 0 and it seems 

that they do not have roots in C - Z, because (27t) ^ 0. 

b) . This forms Gamma(z) = 0 or Gamma(l-z) = 0 do 
not have roots in C-Z. 

c) . Solution of Sin(7r/2z)=0 or Cos(7r/2z)=0. More 
specifically if z = x+yi, then... 

1. Sin(7r/2z ) = 0 Cosh[(^ • y)/2] Sin[(^ • x)/2] + 

+1 • Cos \(tt • x)/2] Sinh[( n • y)/2] = 0 

If k,m g N(Integers) then all the solutions can be 

found with a program by the language of Mathematica. 
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i) .{x = (2+4-k), y = 2-i-(l + 2-m} 

ii) .{x = 4-k,y -2-i-(l + 2-m] 

iii) . {x - (-1 + 4 • k), y = i ■ (-1+4 • m} 

iv) .{x= (-1+4 -k),y = i-(\ + 4-m} 

v) .{* = (1 + 4-&), y = /-(-1 + 4-m} 

vi) . {x = (1+4 • k), y = i ■ (1 + 4 • m} 

vii) .{x = (2 + 4 • k), y = i ■ (4 • m) 

viii) .{x = (2 + 4-k),y - 2-i-(l + 2-m} 

From the generalized solution, it seems that in the pairs 
(x, y) will always arise integers which make impossible 
the case 0 < x <1, therefore there are no roots of the 
equation Sin(;rz/2)=0 in C-Z. 

Also, because as we see, all the roots given by the 
union of the sets x > 1 v x = — 1 v x = 0 all are Integers. 
So it cannot be true that for our case: 

Cos(7r/2z) = 0 Cosh[(yr y)/2] Cos[(/r x)/2] - 

-1 Sin [(n • x)/2] Sinh[( n y)/2] = 0 

If k,m e N{Integers) and the solutions are: 

i) .{x = 4 -k,y = i-(-l + 4-m )} 

ii) .{x = 4-k, y = /-(1 + 4-ra} 

Hi), {x = (-1+4 -k),y = i-(4-m } 

/v).{x = (-1+4- k),y = 2-i-(l + 2-m) 

v) .{x = (l+4-A;), y =4 -i-m} 

vi) . {x = (1+4 • k), y = i - (1 + 4• m] 

vii) .{x = (2+4-k),y = i-(-l + 4-m } 

viii) . {x = (2 + 4-k),y = i-(\+4-m) 

As we can see, again from the generalized solution, it 
seems that in the pairs (x, y) will always arise integers, 
which make impossible the case is 0 < x <1, therefore 
there are not roots receivers of the equation 
Cos(7Tz/2) = 0. Since all the roots are given by the union 
of the sets I>lvi = - lvx = 0it follows that they 
are Integers. 

Great Result. 

«With this three - cases analysis, we have proved that 
the real part of x, of the complex z = x+yi cannot be in 
the interval 0<x<l, and in particular in cases {a, b, c}, 


when they are zeroed. Therefore they cannot represent 
roots in the critical line». 

II) Therefore, now we will analyze the two equations of 
the Riemann zeta function and we will try to find any 
common solutions. 

1. For the first equation and for real values with 
Re(z) > 0 we apply: 

Where f( z ) = 2(27ry z Cos(tt • z/2)r(z) ^ 

=> <^(1 — z) = f(z) 'C(z) but this means that the 
following two cases occur: 

a) . 5(l-z) = <;(z) , where z is complex number. This 
assumption implies that 

C(z) -(l-f(z))=0^> G(z) = o. In theorem 1(1, II), 
(pages 1-2), we showed that the functions £(z) and ^(1-z) 

are 1-1 and therefore if £(x 0 + y 0 i) = g(x' 0 + y'i) 
then , but we also apply that 

C ( x o + 3V) = £> ( x o — )V) = 0 , because we apply it 
for complex roots. The form £(1 — z) = £(z) means 
that if z = X o + y 0 i we apply 1 — x 0 = x 0 namely 

x 0 =1/2 because in this case it will be verified that 

ca -x 0 - y' 0 i ) = ai/2 - y' 0 i) = cat 2 + y’oi) = 

^(1 — x 0 + y' 0 i ) = 0 and because y 0 = y q 

^(v 0 - y 0 i) = £<1/2- y 0 i ) = ^(1/ 2 + y 0 i) = 

= C, (x 0 + y 0 i) = 0 and these are the two forms of the £ 

equation, they have common roots and therefore it can be 
verified by the definition of any complex equation, when 

it is equal to zero. Therefore if z — x 0 + y 0 i then 
x 0 =1/2 which verifies the equation £(z) = 0. 

b) . ^(l-z) + ^(z) , when z is a complex number. 

To verify this case must be: 

C(z) = C( 1 - z) • f(z) =>£(z)= 0 A f(z) = 0 . But 
this case is not possible, because as we have shown in 




Section I(a, b, c), the individual functions of f(z) cannot 
be zero when z is a complex number. 

2. For the second equation and for any real values 
with Re(z) < 1 we apply 

<Z(z) / -z) = 2(2 Try- 1 Sin(7T ■ z i 2)r(l - z) => 
=>az) = Az)-C<X-z), where 


/(z) = 2(2 Try- 1 Sin{n ■ z / 2)r(l - z) but this also 
means that two cases occur: 


a). 5(l-z) = 5(z) , when z is a complex number. 
This case is equivalent to l.a, and therefore if 
Z — x 0 + y 0 l then Xq =1/2 in order to verify the 


# 4.1: 1 st type roots of the Riemann zeta functions (1 st 
equation from the Eq. set). 

For the first category roots and by taking the logarithm of 
two sides of the equations, and thus we get... [5] 

C(1 - z)/f(z) = 2(2 nT z Cos(n • z/2)T(z) => 

=> Logia 1 - z)/C(z)\ = Log\ 2] - zLog\2n\ + 

+ Log\Cos(n • z/2)T(z)] + 2 km 

and the total form from the theory of Lagrange inversion 
theorem, [5] for the root is p x (z) = s which means that 
f(s ) = pi 1 (s) = z> but with an initial value for s which is 

Log [2] 2kri 

Sj ^ --- + - ’ 

Log[2 7T] Log(2 n) 

by setting values for k we can simply calculate the roots 
of the (1 st equation from the Eq. set,(Table.3))... 


equation £(z) = 0 we follow exactly the same process 
algebraically. 


b). 5(l-z)*5(z) , where z is a complex number. 


Similarly, the above case is equivalent to l.b and 
therefore it cannot be happening, as it has been proved. 


#4.Forms of the Riemann £ Functional Equations 

The Riemann £ function has three types of zeros: [7] 

l'C(l - s) = 2 • (2 • n)~ s ■ costa • s/2) • T(s) • ^(s),Res > 0 
Ws) = 2 - (2 tt) s_1 sinta • s/2) • T(1 - s) • C( 1 - s), Res < 1 


Eq. Set And 



(III) 


Usually referred to as the trivial zeros, and non-trivial 
complex zeros.Therefore has been proved that any 
non-trivial zero lies in the open strip 
{s g C : 0 < Re(v) < 1} that is called the critical 
Riemann Strip: And all complex zeros of the function £ 
lie in the line {s e C : Re(v) =1/2} which is called 
the critical line. 

Solving of the Riemann £ Functional -Equations. 

To find the imaginary part we must solve the functional 
equations (Eq. set), and for the cases which the real part 
of the roots lies on the critical line. 

Cases # 4.1 &# 4.2 [6]. 


Therefore, for the first six roots for z we take: 

1st 0.377145562795 

2nd 0.377145562795 +/- 3.41871903296 I 

3rd 0.3771455627955 +/- 6.83743806592 I 

4rth 0.3771455627955 +/- 10.2561570988 I 

5th 0.3771455627955 +/- 13.6748761318 I 

6th 0.3771455627955 +/- 17.0935951646 I 

Table.3-Solution of 1 st equation from the Eq. set, The 

firsts 6 Roots. 

A simple program in mathematica is.. 

k>4;t = (Log[2\ + 2k7ri) I Log[27r]; 

oo 

z = s + ^ ((-1 / Log[27r]) w / Gamma[w +1]) + D[(s') • (. Log[Cos[n • s / 2]] + 

w=1 

+Log[Gamma[s ]] - Log[Zeta[ 1 - s] / Zeta[s]]) w ,{s, w - 1}]; 
N[z/.s-^t,l 0] _ 

With k = 0,±1,±2,... 

But because the infinite sum approaching zero, 
theoretically x gets the initial value 

Log[2] 2*k-7r*i 
Log[2^~] Log(2/z~) 

So we have in this case, in part, the consecutive intervals 
with k = n and k = n+1 for any n >= 4 and for the 
imaginary roots. 

#4.2:2 nd type of roots of the Riemann zeta functional 
equations (2 nd equation from the Eq. set). 

Same as in the first category roots by taking the 
logarithm of two sides of the equations, and thus we get: 
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C(z)/C(l-z) = 2(271 -r 1 Sin(7T ■ Z / 2)r(l - z) => 
Log[£(z)/£(\ - z)] = Log[ 2] + ( Z - 1)Log[2/r] + 

+ Log[Sin(7r • z/2)r(l -z)\ + 2-k-7r-i 
Now, we will have for total roots of 3 groups fields (but I 
have interest for the first group), and therefore for our 
case we will get p x {z) = s which means that: 


1 

2 • n • k , 

Re(z) = — a Im(z) = 

,keZ 

2 

Log[x] 


Following to complete the roots of the two sets of the 
functional equations Eq. Set, we solve the functions as 
cosine or sin according to its Generalized theorem of 
Lagrange [5] 


f(s) = p; l (s) = z, 

but with an initial value for s that is 
Log[;r] 2 • k • n • i 
Sin Log[2;r] + Log(2;r) 

the overall form from the Lagrange inverse theory 
succeeds after replacing the above initial value and 
therefore for the first six roots after calculating them 
we quote the following (Table.4)... 


1st 0.622854437204 

2nd 0.622854437204 -/+ 3.41871903296 I 

3rd 0.622854437204 -/+ 6.83743806592 I 

4rth 0.622854437204 -/+ 10.2561570988 I 

5th 0.622854437204 -/+ 13.6748761318 I 

6th 0.622854437204 -/+ 17.0935951648 I 

Table.4-The firsts 6 Roots of 2 nd equation from the Eq. set 


#4.3: Transcendental equations for zeros of the 
function (Explicit formula) 

The main new results presented in the next few sections 
are transcendental equations satisfied by individual zeros 
of some L-functions. For simplicity, we first consider the 
Riemann-function, which is the simplest Dirichlet 
L-function. 

* Asymptotic equation satisfied by the n-th zero on the 
critical line. [9, 10, 11] 

As above, let us define the function 

x(z) = *~*Pr(z/2 )£(*). 

which satisfies the functional equation 

|x(*) = x(i-*) 

Now consider the Stirling's approximation 


A simple program in mathematica is.. 

k >4;t = {Log[n] + 2kni) / Log[2n]\ 

oo 

z = * + £((-l/ Log[2n]) w / Gamma[w +1]) + D[(s ') • (Log[Sin[n-si 2]] + 

W=1 

+Log[Gamma[ 1 - 5 ]] - Log[Zeta[s ] / Zeta[ 1 - v]])^, {v, w - 1}]; 

N[z!. s-»f,10] 


r (z) = 1 f 2 e~* (1 +0 (z- 1 )) 

where z = x + iy, which is valid for large y. Under this 
condition, we also have 

= exp, (j (y L jg y+ + .rlog y - -y + x + O (y -1 )^ . 


Therefore, with k = 0,±1,±2,... 

But because the infinite sum approaches zero, 
theoretically s gets initial value 

Log[;r] 2 • k • n • i 

z^--—- +- 

Log [2^~] Log (2^~) 

So we have in this case, in part, the consecutive intervals 
with k = n and k = n+1 

for any n >= 4 and for the imaginary roots. And for the 
cases we have for Im(s) the relationship 


2-TT-k 

Im(z) =-, k e Z 

Log[27T] 

For the 3 rd functional equation, one has been previously 
put on the (as in the other two in #4, page 7) and we take 
as imaginary part 


Therefore, using the polar representation 

c = 

and the above expansions, we can write as 
A(x,y) = + i*)| (l + O(^ 1 )j, 

$( x , y) = \ lo S (g^) + - 1) + + iy) + o (y- 1 ) . 


The final transactions we end up with 














10 


Establishing the convention that zeros are labeled by 

positive integers, z n = 1/ 2 + i • y n 

where n = 1, 2,3,4..., we must replace n-^n — 2. 
Therefore, the imaginary parts of these zeros satisfy the 
transcendental equation 



m 

\ + lim — £trg C (| + 5 + -*M = n — If 

<S->0+ 7V 8 

eq. A 

Let us recall the definition used in, namely: 

Si u) = lim — axe 

S-hO-HTT 

:C (j + + iy) = hm [log C (| + 5 + iy)]). 


These points are easy to find, since they do not depend on 
the uctuating S(y). 

We have: 

C(| - m) = C (I + iy) G( y), G{y) - 
The Riemann-Siegel # function is defined by 

&{y) = iug r (-J + fy) - ylog v^r ; 


Since the real and imaginary parts are not both zero, at 
y(+)n then G = 1,whereas at y(-)n then G = -1. Thus 


Q [c(l+^ ) )] =° 

for i?(yi +) ) = (n - 1 )jt, 

^{5+^)] =° 

fonJ(y^ ! ) = (n-g)x. 

they can be written in the form of Eq.b 

f+i 2 jt (ti - 7/8) 

r _ 1 2ir (n — 3/8) 

Vn ^[e-Hn-V 8 )]’ 

~ W [e-!(« - 3/8)]' 


Eq.b 


-CO«f») 



Fig,5 The y(+/-)n are actually the Gram points. 


With this method, we are able to create intervals so as to 
approximate the correct values of imaginary part of 
non-trivial zeros. We did something similar in the cases 
#4.2. What it has left, is to find a method that 
approximates the values of imaginary part accurately. 



The black dots represent the zeros of ^(s) function 
including possible zeros which do not he on the critical 
line (Fig.6). 

The zeros of the Riemann Zeta function. 
“PROGRAMMING” 


#5.1The M Function - Bisection Method... 

Knowing the time of the successive steps (k, k+1) of the 
relationship of imaginary parts 


Im(z) = 


2 -7T-k 

Log[27T ] 


where above ,n = 1,2... and the Labert - Function W 
denotes the principal branch W(0). The y(+)n are actually 
the Gram points. From the previous relation, we can see 
that these points (Fig. 5) are ordered in a regular 
manner... [8,12] 


yl +) < sj ) < d +) < vi 1 < !^ +1 < yi 5 < ■ 


| m ((Eq I, II, III, page 8) with k e. N, we 
Log[7T] 

can calculate the roots by solving the equation 
^(1 / 2 + i • y ) = 0 using the Bisection Method. 
Bisection is the division of a given curve, figure, or 
interval into two equal parts (halves). A simple bisection 
procedure for iteratively converging on a solution, which 
is known to he inside some interval [a,b] proceeds by 
evaluating the function in question at the midpoint of the 
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original interval y = — + ^ and testing to see in which 

2 

of the subintervals [a, (a + b) / 2] or [( a + b) / 2], b] the 
solution lies. The procedure is then repeated with the 
new interval as often as needed to locate the solution to 

the desired accuracy. Let a n , b n be the endpoints at 

the n th iteration (with a x = a and b x =b) and let 

r n be the n th approximate solution. Then, the number 

of iterations required to obtain an error smaller than s is 
found by noting that... h ~ a = -—— and that r is 

J & U n U n 1 U 

defined by r = — (a + b ) In order for the error to be 

J n ^ ' n n 

smaller than s , then... 

\r n - r\ < ^ (b n -a n ) = 2~ n (b - a) < 8 .Taking the natural 
logarithm of both sides gives: 

-/iln2<ln£ — —a) .Therefore, we have for steps: 



#5.2: M-function of the Bisection Method.. 

We define the functions M + d , M d on 
an interval (a,b) according to the scheme: 

I. M d , with M d _ x the Nearest larger of M d _ m 


Mi i + Mi 

—^- d -^,d>2,d>m>d-2 


,, a + b 

Mi =-= 1 

2 

II.M.,with M, , the Nearest smaller of M, 


d-m 


M d -r+M d 


-,d >2,d > m>d —2,M x 


a + b 


,k = 1 


For calculating the roots on solving the equation 
£*(1/ 2 + i • y) = 0 take the limit r 5 according to the 

scheme: 


r = limM^ = limM „,seN with and M ,, 

‘ 3 / 7/7 d d 

n- qo n-oo 


, 2-n-k 7 2-;r-(L + l) N 

(a =-,/? =-), 

Log [2/r] Log [2;r] 


k 


EiV.k>=3 


for the 


x of £"(1/2 + i • ) = 0. 


#5.3 Program in Mathematica for the Bisection 
method of £/l/2 + /• y) = 0. 


Using the Intervals 


2-7T-k ^ _2 - 7r -(k + 1) 
Log\2jz\ ’ Log[2;r] 


and successive steps, we can compute all the roots of 
£Xl/2 + /-y) = 0. 


We can of course use three types such intervals more 
specifically in generall than 


1 ^• n • k b_2’ft-(k + l) 

Log[ 2k] ’ Log[2^-] 


2-n-k , 2-7r-(k + 1) 

-,o =- 

Log[n ] 


3. Were 


y n ^ = y n 


| + . 2tt (n — 7/8) ^ 2 tt (n - 3/S) 

' = H-’lf-Hii-T/S)]’ Vn ' = W[e-i(»-3/8)] ; 


We always prefer an interval that is shorter, in order to 
locate fewer non trivial zeros. 

The most important is to calculate all the roots in each 
successive interval and therefore only then we will have 
the program for data {example: Integer k = 4, and 


8 • n 10 • n , 

- ,b = -) 

Log [27r] Log[27r] 


and Error approximate tol = 


10 A -6 and Trials n=22}.. 


“Programm for Bisection method” 

A program relevant by dividing intervals... 

Clear[ 

f[x_]:=Zeta [ l/2+x*I]; 
k=Input["Epilogh k"]; 
ia=(2*k)*7i/Log[27i]; 

|b=2*(k+l)*7i/Log[27i]; 
tol=Input[ "Enter tolerance " ]; 
n=Input| "Enter total iteration"]; 

|g=N[(-Log[tol]/Log[10]+4)]; 

If[Arg[f[c]]>0,{Print["No solution exists"]}]; 

Print!" nab c . f(c)."]; 

|Do[{c=N[ (a+b)/2,g] ; 

If[Arg[f[c]]<0,a= c,b=c], Print[PaddedForm[i,10],PaddedForm[N[a],{7,7}], 
PaddedForm[N[b],{7,7}],PaddedForm[N[c],{7,7}], 
PaddedForm[N[f[c]],{7,7}]]If[Abs[a-b]<tol && Abs[N[f[c]]]<tol*1000, 
{Print! "The solution is: ",N[c,g]] Exit[]}]},{i,l,n}]; 

Print["The maximum iteration failed,No solution exists"]; 


belong in the interval (a,b),d e N and also 
















This program gives very good values as an approach to 
the roots we ask if we know the interval. Selecting the 
interval for the case 
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, 2 -7T-k 7 2-7T-(k + l). 

(a = - ,b = -) 

Log\2n\ Log[27r] 


k 


e N .k>=3 the results 


are given below in (Table.5).. 


„ , 11 , 

t n =2 7i(n-j)x n 

this equation can be written a 

, -i, 11 , 

x n e " =e («-—) 

Comparing the previous results, we obtain 


1 

Calculation 1st Zetazero zi = — + 14.1347216... 

2 


# 5.5 Programm by Newton’s method, which finds the 
Zeros of g(l/2 + i • y) = 0 

Using Newton’s method we can reach the roots of the 
equation £(l / 2 + i • y) = 0 at a very good initial value 

from y->t n by the explicit formula. 

It follows a mathematica program for the first 50 roots by 
the Newton’s method (Table.6). This method determines 
and detects the roots at the same time in order to verify 
the relation £ (z) = 0 and always according to the 
relation... 



a 

b 

c 

1 

13.6748761300 

17.0935951600 

13.6748761300 

2 

13.6748761300 

15.3842356500 

15.3842356500 

3 

13.6748761300 

14.5295558900 

14.5295558900 

4 

14.1022160100 

14.5295558900 

14.1022160100 

5 

14.1022160100 

14.3158859500 

14.3158859500 

6 

14.1022160100 

14.2090509800 

14.2090509800 

7 

14.1022160100 

14.1556335000 

14.1556335000 

8 

14.1289247500 

14.1556335000 

14.1289247500 

9 

14.1289247500 

14.1422791200 

14.1422791200 

10 

14.1289247500 

14.1356019400 

14.1356019400 

11 

14.1322633500 

14.1356019400 

14.1322633500 

12 

14.1339326400 

14.1356019400 

14.1339326400 

13 

14.1339326400 

14.1347672900 

14.1347672900 

14 

14.1343499700 

14.1347672900 

14.1343499700 

15 

14.1345586300 

14.1347672900 

14.1345586300 

16 

14.1346629600 

14.1347672900 

14.1346629600 

17 

14.1347151300 

14.1347672900 

14.1347151300 

18 

14.1347151300 

14.1347412100 

14.1347412100 

19 

14.1347151300 

14.1347281700 

14.1347281700 

20 

14.1347216500 

14.1347281700 

14.1347216500 


C =- 


2 n(n - —) 
8 


W[e l in - 

O 


where n = 1,2,3. 


Table.5-The first non-trivial Root, Imaginary Part. 

With final value 14.1347216500 , error near of 10 A -7. 
This value is the approximate root of the 
<^(1/ 2 + i • y) = 0 by the nearest error <10 A -7. 

In the event that we have two or more uncommon roots 
within the interval, we divide similar successive intervals 
in the order of finding of the first root either above or 
below. In such a case we have the k = 13, and at the 
interval (47.8620664,51.280785) the two roots are 
48.0051088 and 49.7738324. 


#5.4 Explicit formula and the Zeros of 

4X1/2 + 1-30 = 0. 

Consider its leading order approximation, or 
equivalently its average since 
< Arg£(l/2 + i- y) >= 0 ..[11, 12] 

Then we have the transcendental equation 


—tog (-) = «- — 

In 2m 8 

Through the transformation 


2 n(n - —) 

8 

W[e~\n-^-)] 

O 


where n = 1,2,3.... and where W is the W-function. 

Table[FindRoot[Zeta[s] = 0 , 

{s, 1/2 + N[2*7r*Exp[l] * (n-11/8) /Exp[l] /LambertW[(n- 11/8) /Exp[l]], 
125]*!}] , {n, 50}] _ 


1 

{s->0.5+14,13471) 

11 

{s-> 0.5 + 52.97031) 

21 

(s->0.5+79.3374l) 

31 

(s->0.5+103.726l) 

41 

(s->0.5+124.257l) 

2 

{s-> 0.5 + 21.02201) 

12 

{s-> 0.5 + 56.44621} 

22 

{s->0.5 + 82.9104l} 

32 

(s->0.5+105.447l) 

42 

(s->0.5+127.517l) 

3 

{s-> 0.5 + 25.01091) 

13 

{s-> 0.5 + 59.34701} 

23 

(s->0.5 + 84.7355l) 

33 

{s->0.5+107.169l} 

43 

{s->0.5+129.579l} 

4 

{s-> 0.5 + 30.42491} 

14 

{s-> 0.5 + 60.83181} 

24 

{s->0.5 + 87.42531} 

34 

(s->0.5+1118751} 

44 

(s-> 0.5+131.0881} 

5 

{s-> 0.5 + 32.93511) 

15 

(s-> 0.5 + 65.11251} 

25 

(s->0.5 + 88.8091l) 

35 

)s->0.5+111.875l) 

45 

{s-> 0.5+133.4981) 

6 

{s-> 0.5 + 37.58621) 

16 

{s-> 0.5 + 67.07981} 

26 

(s->0.5 + 92.4919l) 

36 

ls->0.5+114.320l} 

46 

{s->0.5+134.757l) 

7 

Is->0.5 + 40.91871) 

17 

{s-> 0.5 + 69.54641} 

27 

{s->0.5+94.65131} 

37 

(s->0.5+116.227l) 

47 

(s->0.5+138.U6l) 

8 

{s-> 0.5 + 43.32711) 

18 

{s-> 0.5 + 72.06721} 

28 

(s->0.5+95.8706l) 

38 

{s->0.5+118.791l} 

48 

{s->0.5+139.736l} 

9 

{s-> 0.5 + 48.00521} 

19 

{s-> 0.5 + 75.70471} 

29 

(s->0.5+98.8312l) 

39 

(s->0.5+121.370l) 

49 

(s->0.5+141.124l) 

10 

{s-> 0.5 + 49.77381} 

20 

{s-> 0.5 + 77.14481} 

30 

(s->0.5+101.318l} 

40 

(s->0.5+122.947l) 

50 

(s->0.5+143.U2l) 


Table.6- The firsts 50 complex Roots(with W-function) 


A very fast method that ends in the root very quickly and 
very close to the roots of £(s) = 0 as shown.. 
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#5.6 Directly (from Explicit form) with the solution of 
this equation [8] 



Using the initial value the relation 


2n{n - —) 

8 

W[e~\n-—)] 


the explicit formula.. 


As already discussed, the function Argyll 2 + S+i-y) 


according to the intervals defined by the solution of the 3 
based functional equations {#4-1,11,111 p.5}. 

The Riemann Zeta function as defined by Riemann is 
£(s,0)=£(s) and applies the functional equations of 
Riemann.If the equations are solved the only solution for 
the real part of s is 1/2. Many other scientists refer in 
other forms as £(s,a) or £(q(x),a) and many others cases, 
that are unrelated to the original definition and therefore 
relate to the analysis of other cases irrelevant to our 
subject.Is another Hypothesis and not of Riemann. 


oscillates around zero. At a zero it can be well-defined 
by the limit, which is generally not zero. For example, 
for the first Riemann zero yl = 14.1347..the limit 5->0+ 
has value as... 

lim Arg £(1 / 2 + S + / • y) = 0.157873.. 

The arg^ term plays an important role and indeed 
improves the estimate of the n-th zero. 

We can calculate by Newton’s method, and we locate the 
first 30 imaginary part roots, of the equation 
4X1/ 2 + i • y) = 0 ,on the bottom of the (Table.7), 
...x=Im(l/2+yi) 

Table[ 

FindRootfx/(2 tt) t Log[x/(2 tt*«) ] ♦ 1/tt*A rg[Zeta[l/2♦ x*I]] =n-ll/8, 
{x, 1/2 ♦N[2*7T*Exp[l] * (n-11/8)/Exp[l] / LambertW ( (n-11/8) /«) , 20]}), 
{n, 1, 50}] 


1 

{x ->14.1347} 

11 

{x ->52.9703} 

21 

{x-> 79.3374} 

31 

{x -> 103.726} 

41 

{x ->124.257} 

2 

{x ->21.0220} 

12 

{x ->56.4462} 

22 

{x ->82.9104} 

32 

{x ->105.447} 

42 

{x ->127.517} 

3 

{x-> 25.0109} 

13 

{x ->59.3470} 

23 

{x-> 84.7355} 

33 

{x -> 107.169} 

43 

{x -> 129.579} 

4 

{x ->30.4249} 

14 

{x ->60.8318} 

24 

{x-> 87.4253} 

34 

{x -> 111.030} 

44 

{x -> 131.088} 

5 

{x ->32.9351} 

15 

{x ->65.1125} 

25 

{x ->88.8091} 

35 

{x ->111.875} 

45 

{x ->133.498} 

6 

{x-> 37.5862} 

16 

{x -> 67.0798} 

26 

{x ->92.4919} 

36 

{x ->114.320} 

46 

{x-> 134.757} 

7 

{x ->40.9187} 

17 

{x ->69.5464} 

27 

{x ->94.6513} 

37 

{x ->116.227} 

47 

{x ->138.116} 

8 

{x ->43.3271} 

18 

{x ->72.0672} 

28 

{x ->95.8706} 

38 

{x ->118.791} 

48 

{x ->139.736} 

9 

{x->48.0052} 

19 

{x -> 75.7048} 

29 

{x ->98.8312} 

39 

{x -> 121.370} 

49 

{x ->141.124} 

10 

{x ->49.7738} 

20 

{x ->77.1449} 

30 

{x ->101.318} 

40 

{x ->122.947} 

50 

{x ->143.112} 


Table.7-The firsts 50 roots of ZetaZero.(Arg() Method) 


Epilogue ... 

This analysis has proved the Riemann Hypothesis is 
correct, since the real part of the non-trivial zero-equation 
functions is always constant and equals 1/2. This proven 
by 3 independent methods. Theorems {Th l.I-II-pages 
1-5, Theorem 2-pages 5-6,Theorem 3-pages 5-7}. This 
analysis has also demonstrated that the imaginary part of 
the non-trivial zero zeta function accepts certain values 


Part II. Distribution of primes numbers and Zeros 
#6.Indrotuction 

The distribution of prime numbers is most simply 
expressed as the (discontinuous) step function n{x) , 
where 7t(x) (Fig.7), is the number of primes less than or 
equal to v. 



Fig.7 -The primes of function tl(x) 

The function 7i(x) in relation to the random int. x. It turns 
out that 7 r(x) can be expressed exactly as the limit of a 
sequence of smooth functions R n (v). To define R n {x) we 
first introduce the logarithmic integral function Li(x), 



which appears throughout the analytical theory of the 
prime distribution: 

This is a smooth function which simply gives the area 
under the curve of the function 1/log u in the interval 
[2,x]. Don Zagier explains the reasoning behind the 
function Li in his excellent introductory article 
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"The first 50 million prime numbers" [14] , based on his 
inaugural lecture held at Bonn University, May 5, 1975): 
"A good approximation to n(x) , which was first given 
by Gauss is obtained by taking as starting point the 
empirical fact that the frequency of prime numbers near 
a very large number v is almost exactly 1/log v. From 
this, the number of prime numbers up to v should be 
approximately given by the logarithmic sum ... 

Li(x) = 1/log 2 + 1/log 3 + " + 1/log x 

or, what is essentially the same, by the logarithmic 

integral. 



Using Li(v) we then define another smooth function, 
R(x) , first introduced by Riemann in his original 
eight-page paper , and given by 

R{x) = f; ^ Li(x l/n ) 

Riemann's research on prime numbers suggests that the 
probability for a large number v to be prime should be 
even closer to 1/log v if one counted not only the prime 
numbers but also the powers of primes, counting the 
square of a prime as half a prime, the cube of a prime as 
a third, etc. This leads to the approximation: 

tt(j-) + l*(.r ,/2 ) + l?r(x ,/J ) + •••*= Li{x) 

or, equivalently [by means of the Mobius inversion 
formula\ 

tt(x) Li{x) - ~Li{x l/i } - lli(x l/3 }- 


The function on the right side of this formula is denoted 
by R(v), in honour of Riemann. It represents an 


Ii(x ) — 1 T* 


E 


*-i 


(In*)* 
JWt!C(* + l) 


where £(k + 1) is the Riemann zeta 



Fig-8 The non-trivial zeros in the Critical Line 


Here we see the zeros of the Riemann zeta function in 
the complex plane. These fall into two 
categories, trivial and nontrivialzeros (Fig.8).Here are 
some tables on nontrivial zeros compliled by Andrew 
Odlyzko[4].The trivial zeros are simply the negative 
even integers. The nontrivial zeros are known to all lie in 
the critical strip that is 0 < R e[s] < 1, and always come 
in complex conjugate pairs. All known nontrivial zeros 
lie on the critical line Refs] = 1/2. The Riemann 
Hypothesis states that they all lie on this line. The 
difference between the prime counting function and its 
"amazingly good approximation" 

R(x), i.e. the fluctuations in the distribution of primes , 
can be expressedin terms of the entire set of zeros of 
zeta, which we shall represent by p , via the function R 
itself: 

R{x)-n(,x) = Ya R{ - x ^ 

P 


amazingly good approximation to 7l(x) ..For those in Obviously some of the x p are complex values, so here R 


the audience who know a little function theory, perhaps 
I might add that R(v) is an entire function of log x, given 
by the rapidly converging power series: 


is the analytic continuation of the real-valued function R 
defined previously. This was mentioned above by 
Zagier, and is known as the Gram Series expansion: 
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R(x) — 1 -F 


E 


k -1 


(In*)* 

kk\C(k + \) 


show how the jumps in the step function tt(x) can be 
closely approximated by adding to R(x) a correction 
term that involves zeta zeros, we ignore the integral and 


The sum over p separates into two sums, over the 
trivial and nontrivial zeros, respectively. The former sum 
is of course just R(v' 2 ) + R(v" 4 ) + R(v" 6 ) + ', and the latter 
can be written.. 

£[*(**)+ 

k—\ 


the log 2 in second equation; this speeds up the 

computation and will not noticeably affect the graphs, 
especially for x more than about 5.The more zeros we 
use, the closer we can approximate 7t(x) . For larger x , 
the correction term must include more zeros in order to 
accurately approximate n(x) 


The contributions from the complex-conjugate pairs /9 

#7.The GUE hypothesis. 


and p_ K — p K cancel each others' imaginary parts, so 

*(*) = R(x) - £ R(x~ 2m ) + £ T k (x) 

_ m=l _ k= 1 _ 

Where 

r*(*) = -«(*»)-«(. y-‘) 

are real-valued.We can now define the sequence of 
functions R n (x) which approximate n{x) in limit: 


R n (x) = R(x) - E + E T ,*(*) 

-m=l__-*=1-J_ (1 ) 


fix) - 2 Re 


-- 

VEifPilogWJ 4 r—— - di log(2j 

l£r J 1 


■( 2 ) 


th 

where p K is the K complex zero of the zeta function. 

In this formula, E i(z) (Mathematical built-in 
function Explntegral Ei(z )) is the generalization of 
the logarithmic integral to complex numbers. These 
equations come from references [13], [14], and [15]. 
First, let M be the smallest integer such that x llM <2 . 
We need to add only the first M-l terms (that is, 
n = 1,2 ,..M ) in the sum in equation (1). For each of 
these values of N, we use equation (2) to compute the 

value of f(x 1/n ). However, we will add only the first 

N terms (that is, k = 1,2,...N ) in the sum in equation 
(2). Because the purpose of this Demonstration is to 


While many attempts to prove the RH had been made, a 
few amount of work has been devoted to the study of the 
distribution of zeros of the Zeta function. A major step 
has been done toward a detailed study of the distribution 
of zeros of the Zeta function by Hugh Montgomery [18], 
with the Montgomery pair correlation conjecture. 
Expressed in terms of the normalized spacing .. 



this conjecture is that, for M —> oo 

; N + 1 < n < N + M, S n € [a, 0} 


~ f p(0, u)du ~ fl- (1E_^) 2 du . 

Ja Ja mi 

In other words, the density of normalized spacing 

between non-necessarily consecutive zeros is 


1 (Sin miy 
mi 


. It was first noted by the Freeman Dyson, 


a quantum physicist, during a now-legendary short 
teatime exchange with Hugh Montgomery[18], that this 
is precisely the pair correlation function of eigenvalues 
of random hermitian matrices with independent normal 
distribution in (Fig.9&10) of its coefficients. Such 
random hermitian matrices are called the Gauss unitary 
ensemble (GUE). As referred by Odlyzko in [16] for 
example, this motivates the GUE hypothesis which is the 
conjecture that the distribution of the normalized spacing 
between zeros of the Zeta function is asymptotically 
equal to the distribution of the GUE eigen values. Where 
p(0, u) is a certain probability density function, quite 
complicated to obtain for an expression of it). As 
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reported by Odlyzko in [16], we have the Taylor 
expansion around zero.. 

/?(0,w) = — u 2 -2 — u 4 + ... 

3 45 

which under the GUE hypothesis entails that the 
proportion of Sn less than a given small value 8 is 

7T 2 

asymptotic to - S 3 + 0(S 5 ) .Thus very close pair 

9 

of zeros are rare. 


so in particular, for small delta 



so that the probability that the smallest 5n are less than 8 
for M consecutive values of 8n is about.. 

1 “ 0 “ T* 3 ) ^!-exp(-^Af) 

This was the value used in the sixth column of the table. 



Fig9. Probability density of the normalized spacing 

between non-necessarily consecutive zeros and the GUE 
prediction [20] 



The result can be also obtained for the § n + Sn 


7T* 

Prob(<5„ + J n+1 < 6 ) ~ 39 jQQ - JI *» 

from which we deduce the value of the last column 
(Table. 8).. 


Height 

Mini S n 

Maxi S n 

Mini 

+<Wi 

Maxi 
^n+^n + l 

Prob min 
S n in GUE 

Prob min 5 n + 
<$ n+ i in GUE 

10 13 

0.0005330 

4.127 

0.1097 

5.232 

0.28 

0.71 

10 14 

0.0009764 

4.236 

0.1213 

5.349 

0.87 

0.94 

10 15 

0.0005171 

4.154 

0.1003 

5.434 

0.26 

0.46 

10 16 

0.0005202 

4.202 

0.1029 

5.433 

0.27 

0.53 

10 17 

0.0006583 

4.183 

0.0966 

5.395 

0.47 

0.36 

10 18 

0.0004390 

4.194 

0.1080 

5.511 

0.17 

0.67 

10 19 

0.0004969 

4.200 

0.0874 

5.341 

0.24 

0.18 

10 2 ° 

0.0004351 

4.268 

0.1067 

5.717 

0.17 

0.63 

10 21 

0.0004934 

4.316 

0.1019 

5.421 

0.23 

0.50 

10 22 

0.0008161 

4.347 

0.1060 

5.332 

0.70 

0.61 

10 23 

0.0004249 

4.304 

0.1112 

5.478 

0.15 

0.75 

10 24 

0.0002799 

4.158 

0.0877 

5.526 

0.05 

0.19 


Table-8.The values of S n + S n+l 

For very large spacing in the GUE, as reported by 
Odlyzko in [16], des Cloizeaux and Mehta [17] have 
proved that 

logp(0,i) ~ —?r 2 t^/8 (t —* oc) s 
which suggests that 


Fig lO.Probability density of the normalized spacing Sn and the 
GUE prediction [20]. 


max 6 n 

JV+1<ti<JV + jU 


(8 log J/) 1/2 

7T 


#8.Gaps between zeros.. 

The table below lists the minimum and maximal values 
of normalized spacing between zeros Sn and of 

8 n + S n + 1 , and compares this with what is expected 

under the GUE hypothesis. It can be proved that p(0, t) 
have the following Taylor expansion around 0 

P (M = ^ 2 -2^ U 4 +-. 


Statistics of False zetazeros, 8-intervals and 
count of primes. 

#9.General equations. 

With this statistic we find the crowd of individual parts 
intervals defined by consecutive Zetazeros in a fixed 
integer interval. Here we use 8 = 1000. We have the 
general Equation,8 + Zetazero (kin) -Zetazero (kf) = 0 
and after given initial value in kin, i calculate the kf 
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usually by the Newton method. From the initial and 
final value of kin, kf and by performing the process of 
successive intervals, i calculate the number(count) of the 
primes ones that are within the intervals.In this way 
i will have consecutive intervals.. 

S k = {Zetazeroik +1) < p < Zetazero(k),k e ( kin,kf )} 
from where resulting the number of the primes in the 
given consecutive interval. The sum of the primes in 

the given interval 5 = 1000 will be obvious.. 
k f 

s« = U Si 

i—kin 

The False intervals will be in a normalized form 
F s = 1000 - S 6 over N>= 4000. 

#10. The statistics 

The Statistic looking for ways to show us which function 
is the most ideal to get closer the points of interest, uses 
the NonlinearModelFit method of the function 
y = (a + b • x)/(Log(d + c-x) 
and therefore after determining the variables {a, b, c, d} 
we are able to make statistical and probable prediction at 
higher levels of numbers .This function is directly related 
to the function 7t(x) = x / logx which was reportedin the 
introduction that is, defining the number of primes 
numbers relative to x. We will do a double statistic of the 
intervals 

{5, False intervals} and (5, number of primes} that we 
are ultimately interested in the statistics mainly the count 
of the primes inside at 8-intervals. 

#ll.The first statistic is about count of False intervals 
and the count of the primes that corresponding in them. 
The table below gives it in aggregate until the count of 
primes equal p § =29 corresponding to the count of 
False F § =971 and the Total range Integers for a interval 
D=1000 to 2-10 14 + 1000,(Table.9). 


datal= {{499,162}, {735,129}, {825,121}, {885,114}, {887,113}, 
{892, 108}, {895,105}, {894, 106}, {900, 100}, {903, 97}, 


{897, 103}, 

{891, 99}, 

{904, 96}, 

{898, 102} 

, {908, 92}, 

{902, 98}, 

{912, 88}, 

{902, 98}, 

{908, 92}, 

{902, 98}, 

{906,94}, 

{902, 98}, 

{912, 88}, 

{908, 92}, 

{905, 95}, 

{912, 88}, 

{902, 98}, 

{908, 92}, 

{902, 98}, 

{906, 94}, 

{902, 98}, 

{912, 88}, 

{908, 92}, 

{905, 95}, 

{912, 88}, 

{908,92}, 

{914, 86}, 

{900, 90}, 

{908, 92}, 

{912, 88}, 

{925, 75}, 

{932, 68}, 

{936, 64}, 

{930, 70}, 

{937, 63}, 

{934, 66}, 

{945, 55}, 

{925,75}, 

{940, 60}, 

{939, 61}, 

{930, 70}, 

{938, 62}, 

{946, 54}, 

{949, 51} 

{946, 54}, 

{945, 55}, 

{955, 45}, 

{944, 56}, 

{942, 58}, 

{951,49}, 

{956, 44}, 

{962, 38}, 

{953, 47}, 

{956, 44}, 

{955, 45}, 

{961, 39} 

{953, 47}, 

{960, 40}, 

{955, 45}, 

{966, 34}, 

{963, 37}, 

{961, 39} 

{971, 29}, 

{974, 26}, 

{963, 37}, 

{965, 35}, 

{962, 38}, 

{965, 35} 

{970,30} 

{966, 34}, 

{964, 36}, 

{970, 30}, 

{976, 24}, 

{973, 27} 

{971, 29}}; 


Table.9-The data.l 


The diagram given below (Fig.l 1) shows the 
arrangement of the points at the level (x, y) according to 
the datal[19].. 



Fig 11 .The depiction of the intervals False and the primes ones 
(file:datal) contained in the intervals N+5-N = 1000. 

Even more macroscopically (Fig. 12), the points are 
shown by the line y = 988.709 - 0.988372*x, which, as 
they appear, are stacked on its lower right.. 



Fig.12- Graphic depiction of the line y = 988.709 - 
0.988372*x with the archive (data 1) points. 


It is obvious that for {x = 999, y = 1.32} and for {y = 
1000, x = 0.33}, a value located below the unit and 
means that it can in such a interval, and we are talking 
about this is for high order integers intervals so there is 
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not one prime within the interval 5= 1000, chosen at 
random. 


#12..The second statistic refers to the count of the 
primes ones that located in the intervals [N,N+ 5] with N 
= 1000 up to 2 • 10 14 , range 8 = 1000 and the count of 
the primes within successive ZetaZeros .After we 
found the roots kf of general Equation.. 

5 + Zetazero (kin) -Zetazero (kf) = 0 of given initial 
value kin, usually by the Newton method.The datal9 
(Table 10), that we have met with the above method are.. 


data= {{1000,162}, {2000,129}, {3000,121}, {4000,114}, {5000,113}, {6000,108}, 

{8000,105}, {9000,106}, {10000,106}, {11000,100}, {12000,97} 

{13000,103},{14000,99}, {15000,96}, {16000,102}, {17000,92}, 

{18000,98}, {19000,88}, {20000,98}, {21000,92},{22000,98}, 

{23000,94}, {24000,98}, {25000,88}, {26000,92}, {27000,95}, 

{27000,95}, {28000,88}, {29000,92}, {30000,86}, {40000,90}, 

{50000,92}, {60000,88}, {1001000,75}, {2001000,68}, {3001000,64}, 

{4001000,70}, {5001000,63}, {6001000,66}, {7001000,55}, {8001000,75}, 
{9001000,60}, {10001000,61}, {20001000,70}, {40001000,62}, 

{60001000,54}, {80001000,51}, {100001000,54}, {200001000,55}, 

{600001000,56},{800001000,58}, {1000001000,49}, {10000001000,44}, 
{20000001000,38},{30000001000,49},{10000001000,44},{20000001000,38}, 
{30000001000,47}, {40000001000,44},{60000001000,45},{80000001000,39}, 
{100000001000,47}, {200000001000,40},{400000001000,45},{600000001000,34}, 
{1000000001000,37}, {2000000001000,39},{3000000001000,29}, 
{4000000001000,26},{5000000001000,39}, {4000000001000,26}, 
{5000000001000,37}, {6000000001000,35}, {7000000001000,38}, 
{8000000001000,35}, {9000000001000,30}, {10000000001000,34}, 
{20000000001000,36}, {40000000001000,36}, {70000000001000,30}, 
{90000000001000,24}, {110000000001000,27}, {200000000001000,29}}; 

Table.lO-The data 


Using the NonlinearModelFit[19] method of the 
function y = (a + c * x)/ (Log[d + b * x] and after 
specifying the variables {a, b, c, d} the following result 
the function will be.. 

_ 1108.254246288494-1.116325646187134X10“ 12 x 
^ Log[— 4503.9033617702 3+ 5.42 563 5171403 081x] 

with very good approach and value performance in each 
pairing.The diagram (Fig. 13) of the above equation is 
shown in more detail below with gravity in the latest 
data. 



y - (a + c * x)/(Log[d + b * x] .Thus a test value for x = 
2-10 14 gives us y = 25.6, close to 29 we took with the 
analysis.They follow statistics ANOVA and t-Statistic 
[19,20] and we get the (Tables 11,12)... 


Analysis -ANOVA... 



DF SS MS 

Model 

Error 

Uncorrected Total 

Corrected Total 

4 441724. 110431. 

74 1329.65 17.9683 

78 443054. 

77 71143.8 


Table.ll 


Analysis - T-statistic... 



Estimate 

Standard Error 

t- Statistic 

P-Value 

a 

1108.25 

8.79256 

126.044 

4.12076x10-“ 

b 

-1.11633 x10~ 12 

5.63966 x10‘ 15 

-1.97942 

0.0514887 

c 

5.42564 

0.186181 

29.1417 

2.68171 * lO" 42 

d 

-4503.9 

0.153555 

-29330.8 

3.51117 * 10" 263 


Table.12 


As we can see from the results, we have an important 
Standard Error only for the variable a. The other 
variables are observed to have a low statistical error and 
have good compatibility. By adapting the method as we 
see, we associate two lists of results of the number of 
primes and long intervals, which although 
disproportionately, work together with impeccable and 
good contact. 

Ronald Fisher introduced the term variance and proposed 
its formal analysis in a 1918 article The Correlation 
Between Relatives on the Supposition of Mendelian 

Inheritance . His first application of the analysis of 
variance was published in 1921. Analysis of variance 
became widely known after being included in Fisher's 
1925 book Statistical Methods for Research Workers . 
One of the attributes of ANOVA that ensured its early 
popularity was computational elegance. The structure of 
the additive model allows solution for the additive 
coefficients by simple algebra rather than by matrix 
calculations. In the era of mechanical calculators this 
simplicity was critical. The determination of statistical 
significance also required access to tables of the F 
function which were supplied by early statistics texts. 


Figl3.-The graphic depiction of the function 
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#.13. Standar -Error and Confidence -Interval 

Part-1. 

Where we observe, apart from #11-#12 initial 
measurements, the Standar -Error as well and the 
Confidence Interval is stabilized at good and 
acceptable values..(Table. 12) 


Observed 

Predicted 

Standard Error Confidence Interval 

162 

162.352 

4.17901 


154.025, 170.679' 

129 

126.574 

0.851191 


124.878, 128.27 

121 

118.232 

0.759511 


116.719, 119.745} 

114 

113.637 

0.729151 


112.184, 115.09} 

113 

110.529 

0.711538 


109.112, 111.947} 

108 

108.21 

0.699052 


106.817, 109.602} 

105 

104.861 

0.68137 


103.503, 106.219} 

106 

103.582 

0.67462 


102.237, 104.926} 

106 

102.476 

0.668765 


101.143, 103.808) 

100 

101.504 

0.6636 


100.182, 102.826} 

97 

100.64 

0.658982 


99.3265, 101.953} 

103 

99.862 

0.65481 


98.5573, 101.167} 

99 

99.1565 

0.651008 

■97.8593, 100.454 

96 

98.5114 

0.647517 


97.2212, 99.8016} 

102 

97.9179 

0.644293 


96.6341, 99.2017} 

92 

97.3687 

0.6413 


96.0909, 98.6465} 

98 

96.858 

0.638507 


95.5858, 98.1303] 

88 

96.3812 

0.635891 


95.1142, 97.6482} 

98 

95.9342 

0.633431 


94.6721, 97.1964} 

92 

95.5138 

0.631112 


94.2563, 96.7713} 

98 

95.1171 

0.628918 


93.864, 96.3703) 

94 

94.7418 

0.626837 


93.4928, 95.9908} 

98 

94.3858 

0.624859 


93.1408, 95.6309} 

88 

94.0474 

0.622974 


92.8061, 95.2887} 

92 

93.7249 

0.621176 


92.4872, 94.9627) 

95 

93.4171 

0.619455 


92.1829, 94.6514} 

88 

93.1228 

0.617807 


91.8918, 94.3538} 

92 

92.8408 

0.616226 


91.613, 94.0687} 

86 

92.5703 

0.614706 


91.3455, 93.7951} 

90 

90.3458 

0.602126 


89.1461, 91.5456} 

92 

88.7017 

0.592736 


87.5206, 89.8828} 

88 

87.4066 

0.585287 


86.2403, 88.5728} 

75 

71.4688 

0.490798 


70.4909, 72.4468} 

68 

68.4113 

0.47218 


67.4704, 69.3521} 

64 

66.7409 

0.461946 


65.8205, 67.6614} 

70 

65.6044 

0.454957 


64.6979, 66.5109. 

63 

64.7492 

0.449683 


63.8531, 65.6452) 

66 

64.0667 

0.445467 


63.1791, 64.9544} 

55 

63.5009 

0.441964 


62.6202, 64.3815} 

75 

63.0187 

0.438976 


62.1441, 63.8934} 

60 

62.5995 

0.436374 


61.73, 63.469) 

61 

62.2291 

0.434073 


61.3642, 63.0941} 

70 

59.8979 

0.419533 


59.062, 60.7338} 

62 

57.735 

0.405955 


56.9261, 58.5439) 

54 

56.5407 

0.398419 


55.7468, 57.3346} 

51 

55.7229 

0.393243 


54.9393, 56.5064} 

54 

55.1046 

0.389322 


54.3289, 55.8804? 

55 

53.2687 

0.377631 


52.5163, 54.0212} 

45 

51.5512 

0.366632 


50.8207, 52.2817} 


Part-2. 

Continuity... 


Observed 

Predicted 

Standard Error 

Confidence Interval 

56 

50.5969 

0.360493 

{49.8786, 51.3152} 

58 

49.941 

0.356262 

{49.2311, 50.6509} 

49 

49.4438 

0.353049 

{48.7403, 50.1473} 

44 

44.8373 

0.322996 

{44.1937, 45.4809} 

38 

43.6138 

0.314912 

{42.9863, 44.2413} 

47 

42.9283 

0.310357 

{42.3099, 43.5467} 

44 

42.4548 

0.307197 

{41.8427, 43.0669} 

45 

41.8046 

0.302835 

{41.2012, 42.408} 

39 

41.355 

0.299799 

{40.7576, 41.9523} 

47 

41.0126 

0.297474 

{40.4199, 41.6053} 

40 

39.9828 

0.290385 

{39.4042, 40.5614} 

45 

38.9993 

0.283419 

{38.4345, 39.564} 

34 

38.4427 

0.279364 

{37.886, 38.9993} 

37 

37.7577 

0.274288 

{37.2112, 38.3043} 

39 

36.8486 

0.267924 

{36.3147, 37.3824} 

29 

36.3208 

0.265411 

{35.7919, 36.8496} 

26 

35.9441 

0.26516 

{35.4158, 36.4725} 

37 

35.6487 

0.266654 

{35.1174, 36.1801} 

35 

35.404 

0.269629 

{34.8667, 35.9412} 

38 

35.1937 

0.273911 

{34.648, 35.7395} 

35 

35.0085 

0.279366 

{34.4519, 35.5652} 

30 

34.8423 

0.285876 

{34.2727, 35.4119} 

34 

34.6909 

0.293338 

{34.1064, 35.2754} 

36 

33.6014 

0.402929 

{32.7986, 34.4043} 

36 

32.2196 

0.695498 

{30.8338, 33.6054} 

30 

30.6849 

1.16905 

{28.3555, 33.0143} 

24 

29.7968 

1.48821 

{26.8314, 32.7621} 

27 

28.9648 

1.80737 

{25.3635, 32.5661} 

29 

25.5626 

3.23502 

{19.1167, 32.0085} 


Table.12 

#14. Statistical comparison. 

With statistical comparison it appears that as long as we 
moving to a higher order of integer-size that are within a 
given interval, we have chosen 8= 1000, the count of the 
primes diminishes to disappear or to there are 1-2 primes 

at high order intervals of more than 10 20 . This also 

agrees with of Gram’s law. In particular, the problem of 

distribution of the differences t n+1 — t n (that is of 

difference of ZetaZeros) is considered [21]. If we 
accept Gram's law then the order of this difference does 
not exceed the quantity ... 
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C +1 


2;r 

ln(^) 


■> 0, n —^ oo 


for much larger integers [20,21] then mean and their 
mean value is close to,, 


t 


n +1 


lnln(n) 
In (n) 


■» 0,n 


And as the above analysis of - False intervals- as we 
have shown, it is compatible with this result of law 
Gram’s. 


#.15. From above analysis they arise 3 big 
conclusions: 

1 st . The number(Count) of primes located within 8 
intervals is gradually decreasing with a higher order 
size n of 10" . 

2 nd . The distribution of the number(Count) of the 
primes, within interval 8 follows a Nonlinear 
correlation is expressed by the function 
y = (a + b-x) /(Log(d + c-x) similar to the number 
(Count) of the primes 7r(x) « x / Log(x ), which is 
apply approximating for large numbers. 

3 rd .The measurement(Count) of false intervals and the 
count of the primes follows a linear correlation and it 
increases if the order of magnitude size of the integer 
increases. 
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